from sympy import *
from sympy.plotting import plot3d

a, b = symbols("a,b")
eq_1 = 3 * b + 2 * a - 166
eq_2 = 2 * a - b - 67

s = solve([eq_1, eq_2], [a, b], dict=True)
print(s)

size = 100
plot3d(eq_1, eq_2, (a, -size, size), (b, -size, size))

# size = 5
# eq_3 = sin(a) + cos(b)
# # eq_3 = sin(a)
# plot3d(eq_3, (a, -size, size), (b, -size, size))
